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Abstract 

We prove that majorization relations hold step by step in the 
Quantum Fourier Transformation (QFT) for phase-estimation algo- 
rithms. Our result relies on the fact that states which are mixed by 
Hadamard operators at any stage of the computation only differ by a 
phase. This property is a consequence of the structure of the initial 
state and of the QFT, based on controlled-phase operators and a single 
action of a Hadamard gate per qubit. The detail of our proof shows 
that Hadamard gates sort the probability distribution associated to 
the quantum state, whereas controlled-phase operators carry all the 
entanglement but are immaterial to majorization. We also prove that 
majorization in phase-estimation algorithms follows in a most natural 
way from unitary evolution, unlike its counterpart in Grover's algo- 
rithm. 

PACS numbers: 03.67.-a, 03.67.Lx 
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1 Introduction 



Majorization theory emerges as the natural framework to analyse and quantify 
the measure of disorder for classical probability distributions CO 0- Ma- 
jorization ordering is far more severe than the one proposed by standard Shannon 
entropy. If one probability distribution majorizes another, a set of inequalities 
must hold that constrain the former probabilities with respect to the latter. These 
inequalities entail entropy ordering, but the converse is not necessarily true. Quan- 
tum mechanically, majorization has proven to be at the heart of the solution of 
a large number of quantum information problems. It has been shown that ma- 
jorization plays a fundamental role in topics like ensemble realization, conversion 
of quantum states via local operations and classical communication, and charac- 
terization of positive operator valued measurements 

Recently a Majorization Principle has been formulated and checked in the 
efficient quantum algorithms known so far: the family of Grover algorithms and 
the family of quantum-phase estimation algorithms. More precisely, Grover's j!2j 
and Shor's JH] algorithms operate in such a way that the probability distribution 
associated to the quantum state in the computational basis is majorized step by 
step until it reaches an optimal state. This property has been completely proven in 
the case of Grover's algorithm for search of an item in an unstructured database. 
The algorithm is based on the iterative application of a unitary transformation. 
It is easy to see |Hj that each of these transformations majorizes the probability 
distribution associated to the computational basis provided the initial state is 
symmetric. This step- by-step majorization progresses smoothly until the algorithm 
reaches the solution state after 0(y/~N) operations, where ./V is the number of 
entries in the database. 

The case of phase-estimation algorithms, which include Shor's factorisation 
algorithm, Simon's algorithm, clock synchronisation and Kitaev's algorithm for 
the Abelian stabiliser problem ^3] is more subtle. The key ingredients in all these 
algorithms are the use of the quantum Fourier transformation operator (QFT) and 
the promise of a specific structure of the initial state. In 6 it was checked that 
the canonical form of the QFT majorizes step by step the probability distribution 
attached to the computational basis. Here we provide a complete proof of how 
the notion of majorization formulated in |5J explicitly operates in the special case 
of phase-estimation quantum algorithms, and how this notion is related to the 
efficency shown by these algorithms. 

One of our main purposes is to present an explicit and detailed proof of 
the following proposition: Majorization works step by step in the QFT of phase- 
estimation algorithms. The whole property is based on the following ideas: Hadamard 
operators act majorizing the probability distribution given the symmetry of the 
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quantum state, and such a symmetry is partially preserved under the action of 
both Hadamard and controlled-phase gates [BJ. The mathematical formulation of 
these concepts will lead us to define an "H(i)-pair" as that pair of computational 
states that can be mixed by a Hadamard operator acting over the i-th qubit. Fur- 
thermore, we also work out a property concerning the way in which majorization 
emerges naturally from unitary evolution for this type of quantum algorithms. 

In order to assess the significance of step by step majorization we have ana- 
lyzed a series of further examples that will be presented in a separate publication 
[7j. First, we have analyzed a variant of Berstein-Vazirani algorithm [H] where 
a quantum algorithm is able to obtain the slope of a linear function defined on 
Zjv using only one query to the function. The efficient algorithm solving this 
problem is of interest because no entanglement is ever present along the compu- 
tation, while majorization is verified. Second, we have considered the problem of 
determining the parity of a given function using oracle calls as introduced in ref. 
[Sj. The quantum oracular algorithm proposed showed no speed-up as well as no 
step-by-step majorization. Third, we have studied the set of quantum adiabatic 
algorithms proposed by Farhi et al. The efficiency of these algorithms has been 
studied in detail for Grover's problem showing that non-linear evolutions do lead 
to more efficient algorithms We have indeed checked that the efficient time- 
path is associated to step-by-step majorization while the non-efficient time-path is 
not. Finally, we have studied the recently proposed quantum random walk algo- 
rithm to solve a classicaly hard graph problem This algorithm also shows a 
step-by-step minorization-majorization cycle. These four instances of step-by-step 
majorization will be presented elsewhere but can be considered here as motivation 
for a link between majorization and efficiency. 

We have structured the paper as follows. In Sec. 2 we review some concepts 
about majorization theory and how they relate to quantum algorithms. We de- 
velop in Sec. 3 some properties of majorization in phase-estimation algorithms 
and present in detail the main problem to be solved. In Sec. 4 we produce the 
proof for QFT step-by-step majorization. In Sec. 5 we analyse the way in which 
majorization arises in these algorithms and, finally, in Sec. 6 we collect our con- 
clusions. 

2 Majorization theory and its relation to quan- 
tum algorithms 

Let us review the notion of majorization formulated in jS] for quantum algorithms. 
Let us consider two vectors x, y G R d such that Ylt=i x i = Yli=i Vi = 1> whose 
components represent two different probability distributions . We say that distri- 
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bution y majorizes distribution x (written as x -< y ) if, and only if, there is a set 
of probabilities pj and permutation matrices Pj such that 

^ = ^2Pj P jV- (!) 

j 

This definition gives us the intuitive notion that the x distribution is more dis- 
ordered than y, because the former can be obtained from the latter making a 
probabilistically weighted sum over permutations of y. There is an alternative 
definition of majorization which is often more practical. Consider the components 
of the two vectors sorted in decreasing order, written as (z\, ■ ■ ■ Zd) = We say 
that majorizes if, and only if, the following set of inequalities holds: 

k k 

X><X> k = l...d-\. (2) 

i=l i=l 

A third way of defining majorization involves the use of doubly stochastic matrices. 
A real d x d matrix D = (Dij) is said to be doubly stochastic if its entries are 
non-negative, and each row and column of D sums to 1. We say that y majorizes 
x if, and only if, there exists a doubly stochastic matrix D such that 

x = Dy. (3) 

The three given definitions can be proven to be equivalent 4 . 

Majorization can be related to quantum algorithms in the following way: let 
|i/>m) be the pure state representing the register in a quantum computer at an 
operating stage labelled by m = 0, 1 ... M — 1, where M is the total number of 
steps of the algorithm. We can naturally associate a set of sorted probabilities 
p x , x = 0, 1 ... 2 n — 1 to this quantum state of n qubits in the following way: 
decompose the register state in the computational basis, i.e: 

iv> (m) > = E c - m) i*> > ( 4 ) 

x=0 

where {\x) = \x n -i, x n -2 ■ ■ - Xq)} 2 x S^ denotes the basis states in binary notation, 
and x = Y^j=a x fl 3 ■ The probability distribution associated to this state is 

p» = {pW} p M = | C M|2 = |(x|^ m ))| 2 x = 0,1...2 n -l , (5) 

corresponding to the set of probabilities to find each possible output. A quantum 
algorithm will be said to majorize step by step this probability distribution iff 

p<m)^p<m+l) Vm = l,...,M. (6) 
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In such a case, there will be a neat flow of probability directed to the values 
of highest weight, in a way that the probability distribution will be steeper and 
steeper as the time arrow goes on. 

It is important to note that majorization is attached to a probability distri- 
bution defined on a specific basis. Although majorization is basis dependent, the 
basis where the final measurement closing the algorithm is performed is singled 
out. The computational basis is often the natural measurement basis to analyze 
majorization since it gives the probability distribution associated to an eventual 
measurement. We could stop the computation at some arbitrary time and find, 
if step-by-step majorization is holding, that the probability distribution is orderly 
approaching the final one. There may be instances, like the determination of parity 
[Oj, where the measurement basis is different from the naive computational basis. 
The relevant concept of majorization should then be analyzed in the physical mea- 
surement basis. 



3 Majorization in phase-estimation quantum 
algorithms 

Phase-estimation algorithms, initially introduced by Kitaev [TJI|, form a family of 
efficient quantum algorithms |15j that include, for instance, Shor's factoring al- 
gorithm and discrete logarithms Their relevance is due to the exponential 
gain in computational time over known classical algorithms. The basic problem to 
solve can be stated as follows. Given an unitary operator U and one of its eigen- 
vectors \v), estimate the phase of the corresponding eigenvalue U\v) = e~ 2m ^\v), 
4> £ [0, 1) with n qubits accuracy. An efficient solution was found in ^2] and can 
be summarised in the following series of steps, represented by the quantum circuit 
of Fig. 1: 

(i) Prepare the pure state = 1 00 . . .}\v), where |00. . .) will be called the 
register state and \v) is the source state where we have stored the eigenvector of 
the unitary operator U . 

(ii) Apply Hadamard operators 

U H = -^(a 1 + * 3 ) (7) 

over all the qubits in the register state. 

(hi) Apply bit-wise controlled U 3 gates over the \v) state as shown in the Fig 
1, where each U 3 gate corresponds to the application of j times the proposed U 
gate with j = 0, 1 . . . n — 1. 
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|0> - U H 



|o>- 



|0> - C/h 



|«>- 



t/ 2 



|0) + e- 27r ^ 2n " 1 |l) - 



|0) + e- 27r ^ 2n " 2 |l) 



|0) + e- 2 ^ 2 °|l) - U 



fv- 2 



\v) 



Fig. 1: Quantum circuit corresponding to the phase-estimation algorithm. 
Arrows at the end indicate measurements. 




Fig. 2: Canonical decomposition of the QFT operator. By Uj we denote 
the controlled gate |0)(0| + e 2ni / 2J |1)(1|. 
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(iv) Apply the QFT operator 

2 n -l 

QFT:\ q )^^J2 e2ni<!9 ' /2n W) ( § ) 

q'=0 

over the register state. 

(v) Make a measurement of the register state of the system. This provides with 
high probability the corresponding eigenvalue of U with the required precision. 

Let us now go through the algorithm focusing on how the majorization of the 
computational basis probabilities evolves. The application of the Hadamard gates 
in step (ii) to the initial state produces a lowest element of majorization, 

2 n -l 

|^(n) )=2 -n/2 £ | x) | w) ^ (9) 
x=0 
(ii) _ 

yielding to the probability distribution p x = 2 n \/x. The outcome of the con- 
trolled U J gates in step (Hi) is the product state 

|^(m)) = 2~»/ 2 (JO) + e- 2 ^ 2 ""^|l)) • • • (JO) + e- 2 ^ 2 °^|l)) |u) 

= 2_n/2 Exlo 1 e- 2 ^|x)|w) . (10) 

As the action of these gates only adds local phases in the computational basis, the 
uniform distribution for the probabilities is maintained (p x m ^ = 2~ n Vx). 

Verifying majorization for the global action of the QFT is straightforward. 
After step (iv) the quantum state becomes 

2 n -l 

= 2~ n e- 27Tix ^- y/2n) \y)\v) . (11) 

x,y=0 

We then have the probability distribution 

2 

Vy . (12) 

Majorization between the initial (ii) and final (iv) states is verified [Hj, according 
to the definition in ©. The remaining step (v) corresponds to a measurement 
whose output is controlled with the probability distribution py . 

The quantum speed up in the phase-estimation algorithm is rooted in the 
efficient processing of the QFT. It is then essential to investigate whether the 
majorization property so far observed is also present within the QFT step by step. 



E 

x=0 



-2Trix((j>-y /2 n ) 
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4 Step-by-step majorization of the Quantum 

Fourier Transformation in the phase-estimation 
algorithm 

The mathematical statement about QFT we need to prove reads: 
Theorem 

The QFT majorizes step by step the probability distribution calculated in the 
computational basis as used in the phase- estimation algorithms. 

This theorem is seen to emerge from two facts. It is, first, essential that 
the initial state entering the QFT has a certain symmetry to be discussed. Sec- 
ond, the order of the action of Hadamard and controlled-phase gates maintains 
as much of this symmetry as to be used by the rest of the algorithm. More pre- 
cisely, Hadamard gates take the role of majorizing the probability distribution if 
some relative phases are properly protected. Controlled-phase transformations do 
preserve such a symmetry. 

We divide the proof in three steps: the first one will consist on a majorization 
lemma (here is where majorization enters), the second one will consist on a phase 
preserving lemma, and finally the third one will be the analysis of the controlled- 
phase operators in the QFT. As hinted above, we will observe that the only relevant 
operators for the majorization procedure are the Hadamard gates acting over the 
different qubits, while controlled-phase operators, thought providing entanglement, 
will turn out to be immaterial to majorization. 

4.1 A lemma concerning majorization 

Let us introduce the concept of "H(i)-pair", central to this paper. Consider a 
Hadamard gate Uh,i acting on qubit i of the quantum register. In general, the 
quantum register would correspond to a superposition of states. This superposi- 
tion can be organized in pairs, each pair being characterized by the fact that the 
Hadamard operation on qubit i will mix the two states in the pair. Let us illustrate 
this definition with the example of a general quantum state of two qubits: 

|V>) =a|00)+/?|01)+7|10)+<5|ll> 

= (a|00)+/3|01)) + ( 7 |10)+<5|ll)) 

V V 

H(0)-pair H(0)-pair 

= (a|00)+ 7 |10)) + (/3|01)+5|ll)) . (13) 

V v ' " v ' 

H(l)-pair H(l)-pair 
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The second line corresponds to organizing the state as i?(0)-pairs, because each 
pair shares the first qubit value. The third line, instead, organizes the state on 
i?(l)-pairs, because each pair shares the second qubit value. 
We can now formulate the following lemma: 

Majorization lemma 

Let \ip) be a pure quantum state, written in the computational basis as 



(14) 



x=0 



with the property that the probability amplitudes of the computational H(i)-pairs 
differ only by a phase for a given qubit i. 

Then, the probability distribution resulting from Uh,i\iP) majorizes the one re- 
sulting from \ip). 

Proof 

Consider a n-qubit pure quantum state with the assumed property that 
the probability amplitudes of all the H(i)-pairs differ only by a phase for a given 
qubit i. This state can always be written in the following way: 



+ ■ 



\ij>) = oilO.-.O*...) +a 1 e i5l \0...1 i ...) 
+ a 2 n-i|l . . . 0* . . .) + o 2 n-ie W 2"-i |1 . . . V 



(15) 



This expression makes explicit that the amplitude for every pair of states that can 
be mixed by a Hadamard transformation on the qubit i only differ by a phase. 
The Hadamard gate Uh,% will mix all these pairs. The two states in every pair are 
equal in all their qubits except the i-th one. After the application of the Uh,% we 
have 

U H ,i\if>) = 2" 1/2 (ai(l + e iSl )\0 . . . i . . .) + oi(l - e iSl )\0 . . . V . . .) 
+ • • • + a 2 n-i (1 + e'V-i )|l . . . 0* . . .) + a 2 n-i (1 - e*V-i )|1 . . . P . . .)) . (16) 

We have to find a set of probabilities pj and permutation matrices Pj such that 

/ 



/ l«i| 2 \ 
|ai| 2 



|a2»— 1| 
V|a 2 n_i| 2 / 



|ai| 2 (l + cosi5i) \ 
|ai| 2 (l — cos(5i) 



a2 n -i| 2 (l + cos <5 2 n_! 
y|a2«-i| 2 (l - cos 5 2 n-i)J 



(17) 
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and the unique solution to this probabilistic mixture is 



Pi = P2 



n 



Pi 



V 



P> 



V 



/o i 

1 



V 



1 

1 0/ 



(18) 



The permutation matrix P\ is nothing but the identity matrix and Pi is a per- 
mutation of the probabilities of each pair which has undergone Hadamard mixing. 
This completes out the proof. □ 

Note that we have made use of the majorization's definition given in eq. Q 
in this proof. Consequently, it has not been necessary to introduce any notion 
of sorting of the components of the probability distribution. Such a requirement 
would have involved an explicit knowledge of the factors and 5i for all i, which 
is unknown to us in principle. 

The lemma we have proven states that Hadamard transformation do order 
probability distributions when the input state has a special structure, namely those 
amplitudes to be mixed only differ by a phase. This is the key element pervading 
the whole proof. Hadamard transformations and controlled-phase transformations 
carefully preserve such a structure when needed as we shall now see. 



4.2 A lemma concerning phase preservation 



Phase preserving Lemma 

Let us consider a set of Hadamard gates {Uhj} with j = 0, 1 ... n — 1 such that 
each of them can act only once and a quantum state \tp) such that the probability 
amplitudes of the computational H(j)-pairs differ only by a phase 

Then, after applying any Hadamard Un,i operator from this set, the resultant 
state retains the property that the H(j) -pairs differ only by a phase \/j ^= i. 

This lemma states the fact the QFT works in such a way that states to be mixed 
by Hadamard transformations only differ by a phase all along the computation, 
till the very moment when the Hadamard operator acts on it. In other words, the 
structure of gates is respectful with the relative weights of the H(z)-pairs. 

To prove the above lemma we need to build some intuition. Let us first go into 
an example. 
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We start by introducing a new notation for the phases appearing in the source 
quantum state of eq. (fTU)l to be operated by the QFT operator. Let us define 
f3 x = —2ttx4>. Then 



2 n -l 



=2 -»/2 ^ e if3 *\x) . (19) 



x=0 

■vn— 1 



Note that, because x = Y17=o '> we nave that 

n— 1 n— 1 

^ = J] -2vrx i 2V = ^ , (20) 
i=0 i=0 

where = — 27r2 4 0. As an example of this notation, let us write the state in the 
case of three qubits: 

W) = (|000) +e iQ2 |100) +e iQ1 |010) + e^+^lllO)) 

+^72 (|001) +e iQ2 |101) + e iai |011) + e i(Q2+ai) |lll» e ia ° . (21) 

We have factorised the «o phase in the second line. Alternatively, we can choose 
to factorise a± 

| V> = ^572 (|000) + e iQ2 |100) + e ia °\001) + e^+^lOl)) 

+^2 (|010) +e iQ2 |110) + e ia °|011) + e i(a2+ao) |lll» e iai , (22) 

or a 2 , 

1^) = (|000) +e iQl |010) +e iao |001) +e i ( Ql+Q °)|011)) 

+^ (|100) +e iai |110) + e iao |101) +e*( Q1+a °)|lll)) e iQ2 , (23) 

In total, the initial state for three qubits can be factorised in this three different 
ways. This example shows that there are three different ways of writing the quan- 
tum state by focusing on a particular qubit. This property is easily extrapolated 
to the general case of n-qubits: we can always write the quantum state \ifj^ n ^) in 
n different ways factorising a particular phase in the second line. 



Proof 

In the general case we can factorise the ay phase so that the pure state is 
written as 

H) = H {lU) )= ^(\0...W...) + --- + e i ^<**\l. ..&>.. 



2 n/2 
1 
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[\Q...V ...) + ■■■ +e i ^i° Lk \l...V ...))e ia 3 . (24) 



Then, the action of the Uhj transforms the state as follows 

UhM= ^S(|0...0»-...> + --- + e i ^«*|1...0*...>) 

+§=^(|0...P...) + --- + e^^^|l...P...)). (25) 

It is now clear that this resulting state still preserves the symmetry property 
necessary to apply the phase preserving lemma to the rest of qubits i ^ j. The 
reason is that the effect of the operator has been splitting the quantum state in 
two pieces which individually retain the property that all the H(i)-pairs differ only 
by a phase for i 7^ j. If we now apply another Hadamard operator over a different 
qubit i, each of these two quantum states splits in turn in two pieces 







Uh,' 




b) = 








(l+e^Xl+e^Wln 

2 (n + 2)/2 U U -- 


.0*. 


. w . . 


■> + ■■ 




.0*. 


. CP . 


•)) 


1 (l+e^)(l-e-i) nn 

1 2 ("+2)/2 U U - 


. . r. 


. . 0V 


■>+• 


• • + e*^ fc ^ i J Qfc 1 . 


. . r. 


..0> 


••)) 


(l- e ^)(l + e».)fi n 


..0* . 


..V. 


■>+• 


• • + e*^ fc ^ i >J Q!fc 1 . 


..o* . 




••)) 




..IV 


..V. 


■>+• 




..IV 




••)) 



(where we have assumed i > j). The register now consists of a superposition of 
four quantum states, each made of amplitudes that only differ by a phase. Further 
application of a Hadamard gate over yet a different qubit would split each of 
the four states again in two pieces in a way that the symmetry would again be 
preserved. This splitting takes place each time a particular Hadamard acts. Thus, 
all Hadamard gates operate in turn producing majorization while not spoiling the 
symmetry property needed for the next step. This completes the proof of the 
phase preserving lemma. □ 



4.3 Analysis of the controlled-phase operators 

It is still necessary to verify that the action of controlled-phase gates does not 
interfere with the majorization action carried by the Hadamard gates. Let us 
concentrate on the action of ?7/j, n -i, which is the first Hadamard operator applied 
in the canonical decomposition of the QFT. Originally we had 

|V>> = ^72 (|00...) + --- + e^M-i Qfc |01...)) 

+^72 (|10...) + --- + e iE ^-i afc |ll...))e ia "- 1 , (27) 
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where we have taken out the a n _i phase factor. After the action of Uun—l we get 

u H ,n-M = (1 2 tn+i"7 2 1) •••> + ••• + e* 2 **- 1 ° fc |01 • • •)) 

+ (1 2 7n + l)7 2 1) (|10 ■••)+••• + e^Mn-l^|ll .,.)) = |fl) + \b) . (28) 

We repeat the observation made in the second step that the state resulting from 
the action of Ujjn-1 can be- divided into a sum of two states, which we have called 
| a) and \b). For both of these two states the amplitudes of the H(j)-pairs Vj 7= n—\ 
still differ only by a phase. 

We can now analyse the effect of the controlled-phase operators. Notice, first, 
that any two controlled-phase operators acting over the same qubit commute, no 
matter what the phases and the control qubits are. This is a direct consequence 
of the diagonal form of these type of operators, written in the computational 
basis (see for example [El)- Consequently, we can focus on what happens after 
applying a general controlled-phase operator over a given qubit. For simplicity we 
will assume that it will act over the (n — l)-th qubit so it will be applied over the 
state Uu,n-i\tp) (the following procedure is easily extrapolated to the controlled- 
phase operators acting over the rest of the qubits). If the control qubit is the Z-th 
one (I 7^ n— 1), the operator will only add phases over those computational states 
of (|2*5|) such that both the 1-th and the Z-th qubits are equal to 1, so we see that it 
will only act effectively over part of the \b) state. Let us write \b) factorising the 
Z-th phase 

\b)= ^SvjP (|1 ... tf ...) + ••• + e*SM«.n-i °» |1 . . . tf . . .)) 

+ (1 2 7n+i)7 2 1) (|1 1' ...> + ••• + e^.n-i^i. ..!<...)) e <«. . (29) 

It is now clear that we have only added a global phase in the second piece of \b), 
which can always be absorbed in a redefinition of the phase a/. Hence we see that 
no relevant change is made in the quantum state concerning majorization, because 
the amplitudes of the computational H (j)-mixable states Vj 7^ n — 1 still differ 
only by a single phase. The action of controlled-phase operators only amounts 
to a redefinition of phases, which does not affect the necessary property for the 
majorization lemma to hold. We see that the needed phase redefinition can be 
easily made each time one of these operators acts over a particular qubit. 

4.4 Summary of the proof 

We can now collect the three pieces of our argument: 
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• Each Hadamard operation acting on a suitably symmetric state produces 
majorization. 

• Each time a Hadamard operator is applied, the state splits into two states 
preserving the key property of the majorization lemma with respect to the 
remaining Hadamard operators. 

• Controlled-phase operators only involve phase redefinitions, which are im- 
material to majorization. 

It then follows the property that QFT operator majorizes step by step the 
probability distribution in phase-estimation algorithms when the initial state is a 
set of states differing by phases. □ 

We would like to emphasise some relevant points emerging from our proof. 
Controlled-phase operators play no role on majorization, though they provide en- 
tanglement. Curiously, Hadamard operators act exactly in the complementary 
way, providing majorization without providing entanglement. It is also interesting 
to note that the majorization arrow in the quantum algorithm is based on two 
ingredients. On the one hand we have the special properties of the quantum state, 
and on the other hand we have the structure of the QFT. A QFT acting on an 
arbitrary state would fail to obey majorization. In particular a subsequent appli- 
cation of a QFT on the final state obtained above would operate a minorization 
of probabilities, till reaching the original state. 

It is arguable that the proof we have presented depends on the specific decom- 
position of the QFT in terms of individual gates. The underlying quantum circuit 
is not unique and majorization may not be present in alternative decompositions. 
To clarify this point we have analyzed the decomposition of the QFT proposed in 
Ref. [18] . This circuit works in a different way, implementing first a series of gates 
that do not change at all the probability distribution and then a set of Hadamard 
operations that do verify step-by-step majorization. The result remains true that 
an alternative efficient QFT obeys step-by-step majorization. 

5 Natural majorization and efficient quantum 
algorithms 

We now turn to investigate further the way majorization has emerged in the phase- 
estimation algorithm as compared to majorization in Grover's algorithm. This 
comparison was not performed in We shall see that both algorithms work in 
a rather different manner. Majorization is more "natural" in the former than in 



14 



the latter. This naturalness is attached to the absence or presence of off-diagonal 
contributions appearing along the unitary evolution. 

For a search in an unstructured database of a particular item, the best known 
classical algorithm (which is simply the examination one by one of all the items in 
the database) takes asymptotically 0{2 n ) steps in succeeding (where 2 n = N is the 
number of entries) ^21- However, Grover was able to find a quantum mechanical 
algorithm that implements a quadratic speed-up versus the classical one (that is, 
Grover 's quantum algorithm takes asymptotically 0(2 n / 2 ) steps). We will not 
enter into precise details of the construction of this quantum algorithm, and will 
only make few comments on the way it proceeds. 

Grover's algorithm ^2] can be reduced to a two-dimensional Hilbert space 
spanned by the state we are searching \m) and its orthogonal |m~ L )|19j. The 
unitary evolution of the quantum state is given by the repeated application of a 
kernel which amounts to a rotation 

K=( COS0 Q - Sh \ ) (30) 
ysmt/ cosy J 

where cos 6 = 1 — 2/2 n . Other choices of kernels are also valid but the one here 
presented is the optimal one [20]. The initial state of the computation is an 
equal superposition of all the computational states, written as \ip) = 2~ n / 2 |m) + 

1 /2 i 

(1 — 2~ n ) ' \m ) in this 2-dimensional notation. For a given intermediate com- 
putation step the state (a, (5) will be transformed to {a! , [)'). If we wish to express 
the initial amplitudes in terms of the final ones, we have: 

(a\ _ ( a' cos 9 + 0' sin 6\ , . 

Pj ~ \-a'sm0 + P'cos0j ' [ ' 

We now take the modulus squared of the amplitudes, obtaining: 

|a| 2 = cos 2 |a'| 2 + sin 2 # |/?'| 2 + 2 cos sin TZe(a'*p') 
\P\ 2 = sm 2 9 |a'| 2 + cos 2 |/?'| 2 - 2 cos sin Ke(a'* (?) . (32) 

It is clear that if the interference terms were to vanish then majorization would 
follow in a straightforward way from these two relations. To see this, note that in 
such a case we could obtain the initial probability set from the final one according 
to the definition given in eq. Q of majorization, with probabilities p\ = cos 2 9, 
P2 = sin 2 9 and permutation matrices P\ = I and P2 = o\ (the last one is just 
the permutation matrix of the two elements). Majorization would then follow in a 
natural way from just unitary evolution. But this is not the case in Grover's algo- 
rithm, because interference terms do not vanish. Yet it is proven that majorization 
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in Grover's algorithm exists [BJ, although the way it arises is not so directly related 
to the unitary evolution in the way presented here. 

Let us turn back to the majorization in the phase-estimation algorithm and its 
relation to unitary evolution. We write the initial state to be operated upon by a 
Hadamard gate acting over the j-th qubit as 



= c |0...(F...) + c i |0...F...> 

+ --- + C2»-l_j|l...(P'...)+C2»-l|l...l J '...) , (33) 

where we are focusing on the H (j)-mixable coefficients. Applying the Hadamard 
gate over the j-th. qubit we get 



U H M = 2- 1 /2( C0 + Cj .)| . . . & . . .) + 2" 1 /2( Co _ Cj .)| ...V ...) 
+ ■■■ + 2~ 1 /2( C2n _ 1 _ j + C2 „__ 1 )|i . ..&...) 

+2- 1 / 2 (c 2 n_ 1 _ j - C2»_i)|l. ..V ...) . (34) 

Inverting the relations in this last equation, we can find the initial probability 
amplitudes in terms of the final ones. For a given pair of amplitudes c m —j and c m 
we find 

r . - 0-1/2 (o- 1 / 2 (r ■ + r ) + 2 _1 /2 (r -r )) 

c m = 2-V2 (2-1/2 (Cm _ . + Cm) _ 2 -i/2 (Cm _ . _ Cm) ) j (35) 

and making the square modulus we have 

| c .12 _ 1 f IO-1/2 ( . , n| 2 \ 1 / IO-1/2 / ._ x| 2 \ 
I '-m— j I — 2 VI y^m—j \'^m)\ j * 2 \| V^m— J , -m!\ I 



"I - ((Cm— j "I - Cm) (Cm—j C m )) 




-2^e ((Cm-j + C m )*(Cm-j - C m )) • (36) 

As in the Grover's previous example, we observe that if interference terms 
disappeared majorization would trivially arise from this set of relations. In such 
a case, we would only have to choose the set of probabilities and permutation 
matrices of (|18[) to prove this property. But again we note that in general those 
interference terms do not disappear, so majorization does not usually arise in 
this natural way from the unitary evolution (exactly the same that happened in 
Grover's algorithm). For those terms to disappear, there must exist very specific 
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properties for the coefficient q. It is a straightforward exercise checking that the 
interference vanish if and only if 



= a. 



<"m—j 



m-3<- 



(37) 



where a m -j is real. 

It is a remarkable fact that this is the case of phase-estimation algorithms. 
Recalling the lemmas from the previous sections, the interference terms will van- 
ish also step-by-step, so in the case of phase- estimation algorithms, step-by-step 
majorization arises as a natural consequence of unitary evolution, something that 
does not happen in a general case. We also realize that the form of the quantum 
state to be applied to the QFT in eq. 1)10(1 is precisely the unique possible form for 
step-by-step natural majorization to appear. In a way we can say that previous 
steps in the algorithm prepare the state in this unique form. 

Actually the test for natural majorization can be performed for any unitary 
operation acting over any quantum state and consequently for any quantum algo- 
rithm. To make this observation clear, let us write in matrix notation a unitary 
evolution of a general quantum state in an iV-dimensional Hilbert space: 



(38) 



where a, are the amplitudes for the initial state, cJ i the amplitudes for the final 
state, and U corresponds to the N x N matrix of a general unitary evolution. 
Inverting this relation we can obtain the original amplitudes in terms of the final 
ones. Now, if we write W as 















\a' N J 




\a N J 



( 



(39) 



\um ■ • • u NN/ 

we can again calculate the square modulus of the original amplitudes obtaining 



\a\\ 



uu\ 2 Wi\ 2 "I 1" I w ijv| 



a 



N 



+ interference terms 



|ayvT = I^TViPla'il 2 + • 
We see that if the interference terms vanish, we find 



+ l^jVArPla'/vl 2 + interference terms . 















\\(xn\ 2 J 




KM 2 ) 



(40) 



(41) 
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where D is a doubly stochastic matrix such that Dij = \uij\ 2 . Recalling eq. © we 
conclude that the final probability distribution majorizes the original one, in a way 
that majorization arises in a natural way from the unitary evolution. This result 
is reminiscent of other applications of majorization theory in quantum mechanics. 
For instance, ensemble realization is rooted in doubly stochastic matrices of the 
above form: given a density matrix p and a probability distribution pi, there exist 
normalised quantum states \ipi) such that 



if and only if {pi) = DX p , where X p is the vector of eigenvalues of p and = \v,ij 
being the elements of a unitary matrix [£]. 



Our main results in this paper have been to provide an explicit and detailed proof 
that the Quantum Fourier Transformation implements step by step majorization 
in phase-estimation quantum algorithms, along with a comparison between the 
different way in which the notion of majorization operates for the two families 
of efficient quantum algorithms known so far. Further work on other quantum 
algorithms (Bernstein- Vazirani problem, determinaiton of parity, quantum adia- 
batic evolution and quantum random walks for non-trivial graphs) that will be 
presented elsewhere 7_ reinforces the link between step-by-step majorization and 
efficiency. 

The proof of majorization for QFT is based on introducing the concept of 
Hadamard pair, H(i)-pair, and three lemmas. A H(i)-pair corresponds to any pair 
of states which can be mixed by a Hadamard transformation on the bit i. The 
three lemmas show that: 

• A Hadamard gate acting on the bit i of a given register carries out majoriza- 
tion provided the register is made of H(i)-pairs where the two states in the 
pair only differ by a phase. 

• A Hadamard gate on a bit i preserves the relative phase structure of H(j)- 
pairs for j ^ i. 

• Controlled-phase gates do not affect the phase structure of any H(i)-pair. 

Step- by-step majorization then follows since the Quantum Fourier Transfor- 
mation only uses one Hadamard gate per qubit supplemented with immaterial 




(42) 



6 Conclusions 
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controlled-phase gates. The detailed way in which majorization takes place sug- 
gest introducing the concept of natural majorization, defined as a majorization 
which follows from unitary evolution in the absence of off-diagonal terms. Natural 
majorization controls phase-estimation algorithms but is not present in Grover's 
algorithm. A classification of quantum algorithms according to a naturalness cri- 
terium for majorization can be performed. We note that this classification fur- 
thermore corresponds to the well classification in terms of their efficiency. Natural 
majorization might be a distinct feature of exponential quantum speed up. 
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